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Abstract. In the following, we will discuss weighted estimates for the squares of the 
Riesz transforms Rf, 7?^ on L^{W) where W € C^^"^ is an A2 weight. We will show 
that if the "Heat A2 characteristic" of W is sufficiently close to 1 then there is a dimen- 
sional constant c > such that 



||^?|l2,W < 1 + C^[W]aH - 1, 



for all i ~ I, m. This is accomplished by proving a Littlewood-Paley estimate with the 
use of the Bellman function technique. These results follow the ones obtained in [5], [6], 
where similar questions were considered for the case of scalar Muckenhoupt Ap weights 
w. The use of matrix weights turns out to be much more complicated than the scalar 
weights since the environment that we have to work in is no longer commutative. 



1. INTRODUCTION 



Throughout this article, we wih denote W as a, positive, invertible and L^^^^(M'"), d x d 
matrix with complex entries. For such W we define the heat extension to the upper half 
space M™+^ as the convolution of W with the fundamental solution k in 
equation ut — Au = 0, 



pm+l 



of the heat 



W^{x,t) = kt*W{x) 



1 



W{y) exp 



(47rt)T 

We define the heat A2 characteristic of W to be the supremum 



4t 



dy. 



[W] 



sup 



{W\x,t))mW-^f)Hx,t) 



If this number is finite we say that W G Ag. 

The components of the Riesz transform in M™, Rm are defined in the usual way, 

as the convolution operators 



Rif{x) 



vr 2 



Xn 



Im+l 



f{y) dy, 
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i = 1, . . . , m. It is known that for i = 1, ... m we have ||-Ri||2 = Ij where by ||...||2 we denote 
the operator norm m the Hilbert space L^. 

The main theorem of this paper is the fohowing continuity result about matrix weights 
and the square of the Riesz transforms. 

Theorem 1. There is a dimensional constant q > such that for all matrix weights W 
of heat A2 characteristic sufficiently close to 1 the estimate 



i=l ' 

is satisfied, where {ji,j2, ■ ■ ■ ,jn} is an arbitrary subset of {!,... ,m} and a = {aj-}f^i is 
an arbitrary choice of signs. 

Before we discuss the proof of Theorem [T] let us make some comments. It has been 
proven by the last two authors that the same question of the continuity of the operator 
norm of a C alder on- Zygmund operator T in L^{w) with respect to w, 1 < p < 00, has a 
positive answer for scalar Ap weights w (for the Muckenhoupt Ap classes see [1]). That 
is, if the scalar weight w is "close" to the constant weight 1, then the number is 
close to ||T||p. See [5] and [6]. Notice that in our case the heat A2 characteristic allows 
us to have a measure of determining when a matrix weight W is "close" to the constant 
weight Id and that [/djy^h = 1 since ||A;||i = 1. Hence, Theorem [T] is a matrix analogue 
of the results appearing in ^ and |6j. Unfortunately, we only have the result for the 
square of the Riesz transforms and not for the Riesz transforms themselves. In the scalar 
case everything is much simpler since there is a very nice interplay between the Ap classes 
and the BMO space, which allows us to go almost back and forth between weights and 
functions of bounded mean oscillation. Such a relationship does not exist between matrix 
A2 weights and BMO because there is no "good" analogue of BMO for matrix valued 
functions. That is why we have to invent some "new" path connecting closeness between 
W and Id and between the operator norm of an operator in lP'{W) and iP'. The main 
ingredient of the proof is the following Littlewood-Paley type estimate which is obtained 
by using the Bellman function technique. 

Theorem 2. There is a dimensional constant q > such that for all matrix weights W 
of heat A2 characteristic sufficiently close to 1 the estimate 

holds true for all vector functions f,g ^ C^{W^) with values in C^. 

It is known that in order to obtain an estimate as the one presented in Theorem [2] it 
suffices to find a suitable Bellman function. For this reason we also have the following 
theorem that guarantees that such a function exists. 

Theorem 3. For any < 6 < 1, define the domain = {(X,Y,x,y,r,s) E M+ x x 
X X C^""^' X C'^^^'^ : |(x,e)| < X^se,e)^ ,\{y,e)\ < Y^(re,eyi,\ < \\rh^\\ <l + 6}, 
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for all directions e G C^. Let K he any compact subset of Dg. Then there exists a function 
B = Bs^Ki^:Y,x,y,r,s) which is infinitely differentiahle in a small neighborhood of K. 
Furthermore, for any e > 0, B^ ^ can he chosen in such a way that the following estimates 
hold: 

{Ai) 0<S< (l + e)(l + c\/5)xV2yi/2 

(^2) -d^B>2\{dx,dy% 
where c is a constant that depends on the dimension d. 

This function will be obtained as a byproduct of a result first appeared in [4J, which 
states that a similar statement to the one presented in Theorem [T] holds for a version of 
the Martingale transform. More precisely, let us consider a matrix weight W defined on M 
and for each dyadic subinterval / of the real line the operator < W >/= W . It has a 

complete system of eigenvectors e}, which we can assume to be an orthonormal basis 
of and let us also denote by 

hi{x) = -^{xi+{x) - Xi-{x)), 



the Haar function associated to / (we denote by 1+ the right "son" of / and by /_ the 
left "son"). Our Martingale M^, where a = {cr/} is a collection of signs, acts on vectors 
functions / in the following way 



Mf/(x)= <^1if,hie'})2-hjel 



l<k<d 

In [3j it was proven that if the A2 characteristic of W, 

[W]a2 = sup <W>]< W-^ >] 
I 

is sufficiently close to 1 then there is a dimensional constant q > such that we have the 
estimate 

(3) \\M^h,w <l + CdV[W]A2-l- 

Now we have in our disposal all the ingredients of the proof which we are going to present 
in the next section. 

Let us elaborate on the motivation of studying weighted norm inequalities with matrix 
weights. The motivation of studying estimates of this type comes from stochastic processes 
and operator theory. Let us consider a multivariate random stationary process. For simplic- 
ity we consider the case of discrete time i.e. a sequence of d-tuples x{n) = (xi(n), ...,Xd{n)), 
n G Z, of scalar random variables such that E|xj(n)|2 < +00 and the correlation matrix 



Q{n,k) = {Q{n,k)ij}i<ij<d := {Exi{n)xj{n)}i<ij<d, 

depends only on the difference n — k (we use the symbol E to denote the expectation). 
Without loss of generality we can assume that the process is complex valued. It is well 
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known (see [8j) that there exists a matrix valued non- negative measure M on the unit 
circle T whose Fourier coefficients coincide with the entries of the correlation matrix 

Q{n,k) = M{n,k), 

n,k £ Z and that if the process is completely regular then its spectral measure, M, is 
absolutely continuous with respect to the normalized Lebesgue measure m on the unit 
circle, i.e. dM = Wdm. The past of the process is defined as 

Xn = span{x j{k) : 1 < j < d, k < n} 

and the future as 

A'" = span{x j{k) : I < j < d,k > n}. 

By writing "span" we mean the closed linear span in the complex Hilbert space L^(J7, dP). 
If we consider the mapping 

where {ej}i<j<d is the standard orthonormal basis of C^, then we obtain an isometric 
isomorphism between span{x j{k) : 1 < j < d, G and L^(VF). The past and the future 
of the process are mapped to the subspaces of L^(VF) 

Xn = span{z''C'^ :k<n} 

and 

X" = spanlz^C^ : /fc > n}, 
respectively. In this representation the angle between past and future is nonzero if and only 
if the Riesz projection P-^ is bounded in the weighted space L'^(W). All these applications 
and many more aspects of the theory of weighted inequalities with matrix weights are 
thoroughly discussed in the introduction of [9j and the references therein. 

2. THE PROOFS OF THE THEOREMS 

First we begin with the proof of Theorem[3]that claims the existence of a certain Bellman 
function. 

Proof. Using duality in inequality ([3]) we get the following for all / € L'^{W) and g E 
L^{W-^) and J G P 

' ' k 

< (1 + cy^iWU, - 1) ( l^iWf, f)dx) ' ( j{W-^g, g) dx) ' . 

Let us take a look at the sum over k. For fixed / the vectors {efj^k constitute an orthonormal 
basis of and the first summand of the sum over k is the first coordinate of the vector 
A// : = < / — < f >/_ times the first coordinate of the vector Ajg. The same thing 
happens until the last summand. Here we have to notice the following equality, which can 
then be estimated as 
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^ |(< / >7+,ef) - (< / >/_,ef)||(< g >,^,ef) - (< g >7_,ef)||/| 
k 

= Y,\{^ifA)mi9A)\\i 

k 

>|/||^(A,/,ef)(A,5,ef) 



|/| (A,/,A,5) 



So what we actually have proven is that 



T\ ^ \{< f >i+- < f >i-^<9>i+- <9>i-) 



4|J| 



< (l + cV[W^U,-l)( /^W'/) dxy[l{W-^g,g) dx)' . 
Hence, we can define the Bellman function to be the supremum of the left hand side 



BiX,Y,x, y, r,s) =sui>{^ ^ |( </>/+-</>/_,< 5 >/+-< >/_ ) 

=< f >i, y=<g>j,X= f {Wf, f) dx,Y= [ {W-^g, g) dx, r =< W >j, s =< >j} 

J J J J 



over the domain 



DQ = {{X,Y,x,y,r,s)e&x 



X M+ X 

.1 , si 



X C^""^ X C^""^ 



\{x, e)| < X2 (se, e)2 , e)| < (re, e) 2 , 1 < ||r2s2 1| < 1 + S}, 

for all directions e in and a fixed 6 G (0, 1). 
Obviously, the function B satisfies 

1 , ,1 



(4) 



0< B{X,Y,x,y,r,s) < {l + cV5)X^Y^ 



Notice that B does not depend on the choice of the interval J. This happens because 
averages of functions remain the same under the use of a linear transformation. In addition, 
we claim that for all 6-tuples a"*" = y+, x"*", y"*", r"*", s"*"), a~ = {X~ ,Y~ ,x~ ,y~ ,r~ , s~) G 
Dg, such that € Dg, the following inequality is true 



(5) 



B 



a+ + a-\ B{a+) + B{a-) 1 



> -X ,y^ -y~ 
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To prove this claim let us consider a positive e. There exists functions f'^ , on J_|_ 

such that they satisfy the conditions in the supremum of the function B for the vector a"*" 
and 

B{a+)-e<-^ \{<f+>j,-<f+>j_,<g+>,^-<g+>:_)\\I\. 
' I&V{J+) 

Do the same for the vector a~ in the interval J_. Define the functions F, G, W on the inter- 

valJas:F=|^+ G = ['^ ' °" , and, W = {^^ 

I /_ , on J_ Iff-, on J_ I W- , on J_ 

Observe that they satisfy the required equalities in order to be acceptable for the supremum 

'+c 
2 



that defines the Bellman function for the vector " and therefore, 



1 1 



E \ift-fl,9t-9t)\\I\ 



4 2|J,, 

+ ^ |(/7+-/7_,57+-57Jim + ^term(/= J) 

> i(B(a+) - e) + i(i?(a-) - e) + \\(x+ - - y.)\. 

Now we need to mollify this function B, in order to make it smooth. This can be done 
in exactly the same way as in [2]. The concavity inequality remains the same after the 
mollification and the size condition can become l+C^e times worse, where Ck is a constant 
that depends on the compact set K. □ 

The next step is to prove Theorem [2j But, now that we have the Bellman function this 
proof is basically the same with the one given in f7]. The only difference is that we work 
with the function 

v{x, t) = {{Wf, f)\x, t), {W-'g, g)\x, t),f\x, t),g\x, t), W\x, t), iW-')\x, t)). 

Finally, for the proof of Theorem [1] we need to observe only the following result (for a proof 
see [2J). 

Proposition 4. Let ^ be real functions in C^(R™'). Then the integral J • dxdt 
converges absolutely for all i = 1, . . . , m and 

RU ■ Hx = -2 [ ^ ■ ^ dxdt. 

J dxi dxi 

Therefore, the left hand side of inequality ([T]) is bounded from above by the left hand 
side of inequality ([2]) and this means that we have the desired estimate. 
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